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ABSTRACT 

Several recent studies show that bright, intermediate and high redshift optically 
and radio selected QSOs are positively correlated with nearby galaxies on a range 
of angular scales up to a degree. Obscuration by unevenly distributed Galactic dust 
can be ruled out as the cause, leaving weak statistical lensing as the physical process 
responsible. However the amplitude of correlations on <1° angular scales is at least 
a factor of a few larger than lensing model predictions. A possible way to reconcile 
the observations and theory is to revise the weak lensing formalism. We extend the 
standard lensing formulation to include the next higher order term (second order) in 
the geodesic equation of motion for photons. We derive relevant equations applicable in 
the weak lensing regime, and discuss qualitative properties of the updated formulation. 
We then perform numerical integrations of the revised equation and study the effect 
of the extra term using two different types of cosmic mass density fluctuations. We 
find that nearby large-scale coherent structures increase the amplitude of the predicted 
lensing-induced correlations between QSOs and foreground galaxies by ~ 10% (not a 
factor of several required by observations), while the redshift of the optimal, i.e. 'most 
correlated' structures is moved closer to the observer compared to what is predicted 
using the standard lensing equation. 



Introduction 



Weak gravitational lensing manifests itself in two observable ways: if a source is extended, 
like a galaxy, its image will appear to be sheared tangentially with respect to the foreground mass 
concentration. If sources are unresolved, like QSOs, weak lensing can be detected in a statistical 
way by the sources' angular (anti) correlation with the tracers of the intervening mass distribution. 

Statistical QSO-galaxy associations on large angular scales, 10' — 1° are the subject of the 
present paper. The mechanism responsible for the associations is believed to be the magnification 
bias, which depends on the shape of the number-magnitude counts of the background sources. If 
the number counts are a power law with a slope a = d\ogN{m)/dm then the over- or underdensity 
of sources down to a certain limiting flux behind the lens is related to the lens magnification M by 
q(M,a) = M 2 ' 5a_1 . The factor 2.5a in the exponent accounts for the magnification of individual 
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sources brightened into a flux limited sample, while the factor M _1 corrects for the area dilution of 
the number density of sources on the sky. Thus if a > 0.4 (a < 0.4) correlations (anticorrelations) 
are predicted between galaxies and background sources. 

There are a number of observations of correlations between various types of QSOs and 
foreground galaxies. Here is a sampling of the recent literature. Radio selected 1 Jy QSOs at 
z s > 0.5 (Kiihr et al. 1981) are correlated with the APM Catalog galaxies, zi ~ 0.1 — 0.3, (Seitz 
& Schneider 1995, Bem'tez & Martinez-Gonzalez 1995, 1997, Norman & Williams 1999) on scales 
between a few arcminutes and a degree. Bright optically selected QSO from the LBQS Catalog 
(Hewett et al. 1995) are correlated with the APM galaxies on degree scales (Williams & Irwin 
1998), while faint UVX selected QSO candidates are anti-correlated with galaxy groups on <10' 
scales (Croom & Shanks 1999). The amplitude of (anti)correlations varies between 30% and 
1 — 2% depending on the angular scale and the details of each study. 

As judged by their redshifts, QSOs and galaxies in all these studies are unconnected, therefore 
there are two possible physical reasons why these would appear to be associated in projection. 
One is weak lensing, as we have just described, the other is dust obscuration. However, a number 
of recent studies, if considered together, reject dust as the likely explanation. 

Depending on the limiting magnitude of the QSO sample magnification bias will induce 
foreground galaxies to be either positively or negatively correlated with QSO. Bright end of the 
QSO number counts has a steep slope, a > 0.4, and so q > 1; the opposite is true for QSO samples 
with faint limiting magnitudes, where a < 0.4. And in fact q is observed to be > 1 in most of the 
above listed studies because most of them use bright QSOs. As one goes to fainter QSOs q seems 
to drop to ~ 1, while faint QSOs, those in the Croom Sz Shanks study are anticorrelated with the 
intervening mass, i.e. q < 1. However to explain q > 1 and q < 1 cases with dust obscuration one 
needs to invoke two different types of dust: Galactic dust will obscure both QSOs and galaxies 
in some directions thus leading to positive correlations between the two populations, while dust 
intrinsic to lenses, i.e. galaxy clusters/groups will obscure QSOs in the direction of the lenses only, 
leading to anticorrelations. This double-type effect makes dust a rather contrived scenario. 

Furthermore, if Galactic dust is responsible for the observed positive correlations one would 
expect the radio selected QSO samples to be less affected, and thus show smaller cross-correlation 
signal compared to optically selected QSOs. The observational situation is exactly opposite: 1 
Jy QSOs are stronger correlated with APM galaxies (Norman & Williams 1999) than are LBQS 
QSOs (Williams & Irwin 1998) on similar angular scales. 

Qualitative evidence not only rejects dust as the underlying physical cause, but also points 
to lensing. Amplitude and significance of correlations generally increase with brighter QSOs, as 
predicted by magnification bias. QSOs at intermediate redshifts, z s ~ 1 — 2, i.e. those roughly at 
the optimal distance for being lensed by structures at z\ ~ 0.1 — 0.3 show the strongest correlations. 
Dust obscuration should show no preference for any particular QSO redshift range. 

All evidence combined strongly suggests that lensing is responsible for the observations. 
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In fact, lensing would have been long accepted as the cause were it not for the unexpectedly 
large amplitude of the (anti)correlations. On small angular scales, ^10' the non-linear growth 
of mass fluctuations with cosmic time produces the correct correlation strength; Dolag & 
Bartelmann (1997) and Sanz et al. (1997) reproduce the 1 Jy-APM correlations observed by 
Bem'tez & Martinez-Gonzalez (1995, 1997). However, for correlations on larger scales, <1°, and 
anticorrelations on <10' scales (Croom & Shanks 1999) the strength of associations exceeds 
predictions by about a factor of 5-10. 

There are three avenues within the lensing hypothesis for reconciling the model predicted 
and observed correlation amplitudes. One is to make the slope of the QSO number counts in the 
appropriate redshift range very steep. Williams & Irwin (1998) estimate the required slope a to be 
about 8. Because weak lensing induces small magnifications, the number counts slope of 'lensed' 
QSO, i.e. those seen in the directions of intervening mass concentrations should be almost as steep 
as the slope of the 'unlensed' QSOs. In other words, the overall QSO number counts slope should 
also be around 8. This is grossly inconsistent with the observed slope of ~ 1.1 — 1.6. We consider 
this option the least likely of the three. 

The second alternative is to invoke mass density fluctuations on <2 — 10/i _1 Mpc scales 
(10' — 1° at zi ~ 0.1 — 0.3) that are substantially larger than observations seem to indicate. 
Williams Sz Irwin (1998) and Croom & Shanks (1999) independently estimate that Q,as should be 
about 3-4 to explain their respective results. The commonly accepted value for Q°' 6 as is ~ 0.6, 
and is derived by two different methods: bulk flows in the nearby Universe and the abundance of 
rich galaxy clusters (Branchini et al. 1999, White et al. 1993). While scale-dependent biasing may 
allow one to increase the amplitude of mass fluctuations by 10 — 50%, an increase by a factor of 
6-8 is outside the range of acceptable possibilities. 

The third option is that the weak magnification regime of the lensing theory needs a revision. 
If future surveys support the results of present observations then this option will have to be 
considered seriously. 

In principle many aspects of the currently accepted description of light propagation in an 
inhomogeneous medium can be questioned, for example the applicability of the Friedmann- 
Robertson- Walker metric and its weakly perturbed version in cosmological context, the validity 
of the general relativistic geodesic equation of motion as applied to light rays, etc. But the 
overwhelming success of the standard picture, which incorporates these concepts, in describing our 
Universe in general and observations of most types of gravitational lensing in particular suggests 
that the modifications to the standard description of lensing has to be sought within the existing 
framework. The standard lensing equation is the solution to the first order approximation of the 
full geodesic equation of motion for photons. There have been attempts in the literature to use 
the full geodesic equation to propagate light rays through a clumpy universe (Holz & Wald 1999, 
Tomita et al. 1999 and references therein), however these studies were focused on small angular 
scale lensing effects, and the equation was used in conjunction with a limited range of prescriptions 
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for mass fluctuations, and thus not fully explored. 

Instead of taking on the full geodesic equation with all the higher order terms, our aim in 
this paper is to go one step beyond the standard lensing description, i.e. study the effect of the 
second order term(s) in the geodesic equation. The advantage of this approach is that we can 
derive useful analytic results, for example, an approximate form of the revised lensing equation 
(Section 3), and the equation for the increment in source magnification arising from retaining the 
second order term(s) (Section 4), and discuss qualitative predictions of the revised formulation 
(Section 5). With the help of numerical integrations we examine how the effects of the second 
order term(s) depend on the type of mass fluctuations populating the universe (Section 6 and 7). 



2. Light propagation equation with second order terms included 

The derivation of the standard lensing equation can be found in many places in the literature 
(Schneider et al 1992, and references therein; Dolag k, Bartelmann 1997); here we follow the 
derivation presented by Kaiser (1998), but for simplicity adopt a spatially flat cosmological model. 
Unlike Kaiser and other authors we do not truncate the geodesic equation at the first order terms, 
but retain the second order term(s) as well. 

The FRW spacetime metric with superimposed weak perturbations is given by 

ds 2 = gapd^dr 13 = a 2 (r/)[-(l + 2$)*? 2 + (1 - 2<5>){dx 2 + dy 2 + dz 2 )}, (1) 

where a{rf) is the scale factor, <3? is the Newtonian potential, r] is the conformal time related to 
cosmic time, t though dt = adr], and dx, dy, dz are the comoving separations between two adjacent 
points in space. Photon trajectories are solutions of the geodesic equation, 

d 2 r a aB{ 1 \dr v dr^ 



9 {9^-o 9 ^)^\^\> (2) 



dx2 v 2 *„wj dX dX , 

where A is the affine parameter, and r ' 1,2 ' 3 stand for rj, x, y, and z (redshift is designated by z 
with a subscript, either z\ or z s ). The elements of the diagonal metric tensor g a p can be read off 
from the line element, eq.(l). Let us write down the x-component of the geodesic equation keeping 
terms up to second order in small quantities, i.e. <&, dx/dX, and dy/dX: 

(fx ^ d$/dr]\2 d$/dz\2 d®/dx\/dz\ 

dtf~~dx~\dX~) "^VdXJ + ~d^\dx)\dx)' ^> 

The first two terms in this equation are the standard first order terms; the last term is the only 
surviving second order term. Changing the differentiation throughout to that with respect to z 
and using that drj/dz = 1 — 2$ to first order, we get 
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where the dots represent differentiation with respect to z, and x is the transverse comoving 
separation between a ray and the z-axis. A similar expression can be written for the y-component. 
This is the equation of motion for the x-component of a single light ray as it wonders through a 
spatially flat universe populated with small mass perturbations. The standard equation (see eq.[8] 
of Kaiser 1998) is eq.(4) minus the last term. 

An alternative way to arrive at eq.(4) is by starting from Fermat's Principle which states that 
the path taken between two points in space, P\ and P2, is a local extremum, 

[Pi 

5 dr] = 0. (5) 
JPi 

From the metric of eq.(l) and the condition that light rays are null geodesies we get 

dn = {l-2$)(dx 2 + dy 2 + dz 2 ) 1/2 , (6) 

which can be rewritten as 

n|t| = 1, (7) 

where n = (1 — 2<3?) is the refractive index of the space-time described in eq.(l), t = xx + yy + zz 
is the unit velocity vector of the light ray, and the dot represents d/drj. Equation (5) can now be 
rewritten as 

[Pi 

5 / n\t\drj = 0. (8) 
JPi 

The integrand is an extremum when Euler-Lagrange equations are satisfied, so in the x-direction 
we must have 

d r3(n|t|)i 9(n|t|) 



dx 



dx 



dr\ 

which is simplified to a differential equation for light rays: 

d_ 

drj 

Rewritting this as 



0, (9) 



nt 



- Vn = 0. (10) 



d 

nt = Vn — t—n, (11) 
dr\ 

makes it apparent that we have arrived at the vector form of eq.(4). (This is the same as eq.[4.17] 
of Schneider et al. [1992], without the last term in their equation.) The term Vn is the gradient 
of the refractive index evaluated in the coordinate system whose z— axis is aligned with the 
unperturbed light ray, while the term tdn/dn is along the direction of the actual ray. Thus the 
RHS of eq.(ll) is the gradient of n in the plane perpendicular to the direction of the actual ray. 

Because the last term in eq.(ll) (i.e. the second order term) is small compared to the first 
one, is it neglected in the standard lensing equation. In physical terms this amounts to setting 
x and y to zero in eq.(ll), so that t is along the z— axis. In other words, this approximation, 
sometimes called the Born approximation, assumes that the actual ray is always parallel to the 
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unperturbed ray. In this paper we use the full eq.(ll) and thus test the validity of the Born 
approximation. 

Using (4) we now write down the comoving separation between two adjacent rays, say rays 
on the opposite sides of a small light bundle: 

Ax « -2Ax- T - 2Ay— - + 2— As, (12) 
Ox z oxoy Oz 

and similarly for the y-component. Here, the terms 2Ax-^^x and 2AyJ^x are much smaller 
than the third term on the right hand side in the above expression, so they were omitted. 



3. Approximate solution: revised lensing equation 

Equation (4) without the second order term has an exact solution, 

rz f)<b 

x = 9 z-2 —(z-z')dz', (13) 
Jo ox 



where 6o is the observed angle of the source with respect to the z-axis in the absence of 
perturbations. If all the mass fluctuations between the observer and the source are confined to a 
plane perpendicular to the z-axis then (13) reduces to the standard lensing equation. There is no 
useful solution for the full equation (4). To find an approximate solution in the weak lensing limit 
let us assume that x is constant: 

x « x/z « Oq. (14) 

This is not a bad approximation in the weak lensing case where path deviations suffered by light 
rays are small. With this assumption the second order term in eq.(4) can be integrated by parts, 

Thus an approximate solution to equation (4) is 

r z <9<E> r z <9<l> 

x^e z-2 —(z-z')dz' + 2x —(z-z')dz'. (16) 
Jo ox Jo oz' 

This is the revised version of the '3-dimensional' lensing equation (13). The equation for the 
transverse separation between two adjacent rays is then, 

. „ r n [ z (d 2 <5> d 2 § \ z'{z - z') „ f z d<5> (z - z') , ,i 

Axn A6 z 1-2 + ^ J -dz' + 2 —± dz' , 17 

L Jo ^ox z oxoy) z Jo oz' z J 

where we have approximated Ax by A#o, and set Ax ~ Ay ~ A9oz' , similar to eq.(14), i.e. 
angular separation between two adjacent rays is constant to zeroth order. 
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4. Magnification of sources 

If the source does not suffer much distortion and the x and y-axes are along the major and 
minor axes of the image, then the magnification of the image relative to the source is approximately 



M 



(^)(^)- <*> 



Ax > V Ay / 

The standard magnification in the weak lensing regime is then given by 

r z /d 2 & <9 2 <E>\ z' 

M -«" i+i L^ + w)' (1 -^ )iz <19> 

and the difference between M st( i and magnification with the second order term is, 



<9$ , z 



AM = M 2nd - M std « -4 / _(!-_) dz'. (20) 



dz 

This additional magnification depends on the gradient of the Newtonian potential along the z-axis, 
or along the line of sight to the source. Equation (20) was found to be in excellent agreement 
with the results of numerical integrations (see below) for magnifications M st d<1.3. This range is 
more than sufficient for the present study; magnifications typically produced by smoothed large 
scale structures on <1° scales are less than 1%. For example, mass fluctuations associated with 
the 18.5 < tur < 20.0 APM galaxies in Williams & Irwin study are of the order 1.005. 



5. Properties of the equation 

Qualitative consequences of the second order term in the geodesic equation can be assessed 
from eq.(17) and (20). Let us consider a single positive mass fluctuation between us and the 
source, d^/dz will be negative between us and the lens, and because of the (1 — z' / z) weighting 
in the last term of eq.(17) it will decrease Ax, i.e. increase magnification. The opposite is true 
for fluctuations of negative Sp/p. Alternatively, looking at eq.(ll) we see that the second order 
term increases the bending of the rays around positive and negative potential wells beyond what 
is predicted by the standard first order term. Thus the second order term reinforces the amplitude 
of (de) magnification predicted by the standard formulation. 

According to eq.(17) both Ax and Ay-components of the beam's cross section are affected in 
the same way by the second order term, i.e. a factor proportional to Jq §^ (1 — y) dz' is added to 
both, therefore magnification along the x-direction is not at the expense of magnification in the 
y direction, and so for small magnifications there will be no change in the lensing-induced shape 
distortion of the image. This was verified by the numerical integrations of eq.(12). 

When a significant lens is encountered, meaning the projected mass density fluctuation is 
high, the middle term in eq.(17) becomes large, as it is related to the mass density through the 
Poisson equation. Thus the magnification is dominated by the second derivatives of and the 
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extra magnification supplied by the second order term, AM is negligible compared to {M st d — 1). 
In other words the standard equation is completely adequate in the regions of high density, for 
example in galaxy clusters and near individual galaxies. 

An important aspect of the standard lensing equation is that not all equal amplitude mass 
fluctuations between the observer and the source contribute equally to the source magnification; 
eq.(19) shows that the fluctuations are weighted by z'(l — z'/z), i.e. lenses approximately half way 
between the source and the observer are most influential, d^/dz in eq.(17) and (20) is weighted 
by (1 — z'/z), thus the lenses closer to the observer contribute more through AM. In practice it is 
probably the combination of the amplitude of fluctuations, their physical scale and proximity to 
the observer that determines their relative contributions to AM. For example, because (1 — z'/z) 
depends on the distance only weakly high frequency fluctuations will mutually cancel out, and will 
not contribute significantly to AM. AM is probably dominated by the 'last screen' of large scale 
structure, like superclusters, voids and filaments, located closer to the observer than the standard 
optimal lensing distance. 

This is an interesting point in the light of some of the observational results. A few 
cross-correlation studies seem to indicate that the galaxies most strongly correlated with the 
QSOs are not the ones at the optimal lensing redshifts but are somewhat closer than that to the 
observer. This was stressed by Bartelmann & Schneider (1993) who reanalyzed Fugmann (1990) 
result of correlations between Lick Catalog galaxies and 1 Jy sources. It was also found to be 
the case by Williams &; Irwin (1998) who show that correlations are strongest for zqso > 1- The 
optimal lens location in the Einstein-de Sitter universe model for sources at zqso > 1 would be 
0.4, while the galaxies in their study lie at z\ ~ 0.2. 

Equation (20) for AM shows that flux is conserved with the addition of the second order 
term, as d^/dz will average out to zero for a collection of sources. 

Finally, because the second order term depends on d^/dz it is essential that the continuous, 
3-dimensional nature of the mass fluctuations is taken into account. If the effect of the second order 
term proves to be non-negligible then the treatment of weak QSO lensing based on multiple-plane 
formalism, where the standard lensing equation is solved successively on many 2D planes stacked 
perpendicular to the optical axis, should not be used. 

We will return to some of the qualitative remarks made in this Section in Section 7, where we 
discuss the results of numerical integrations of eq.(12). 



6. Application of the equation 
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6.1. Numerical integrations 

To study the effect of the second order term further we need to consider an inhomogeneous 
mass distribution in the Universe, solve eq.(12) for many lines of sight and compute relevant 
output quantities, like distribution of magnifications, average magnification, etc. Given the 
simplicity of the approximate solution, eq.(17), this can be done analytically for simple mass 
distributions. For an arbitrary mass distribution a numerical integration approach is preferred. 
We use 5th order Runge-Kutta method to integrate eq.(12) and its y-component counterpart. In 
principle, as a light bundle propagates through a clumpy medium the path of its central ray will 
deviate from the fiducial ray in the absence of perturbations; however, one can ignore this effect 
in the weak lensing regime and evaluate the gradients of the Newtonian potential in eq.(12) along 
the fiducial central ray. When integrating eq.(12) we use boundary conditions at the observer; for 
every line of sight we start with two images seen by the observer, each one confined to the x- and 
y-axis entirely, and propagate these backwards to the redshift of the source. The shape, size and 
orientation of the two final sources per line of sight allows us to solve for the four components of 
the total magnification matrix. All sources are located at z s = 1.5, typical redshift of QSOs in 
correlation studies. 

6.2. Mass density fluctuations 

We are primarily interested in how the weak lensing model predictions are affected by the 
inclusion of the second order term, in particular if the differences in the standard vs. revised 
predictions depend on the type of mass density fluctuations present in the universe. Therefore we 
pick two rather different, but simple mass fluctuation scenarios: randomly distributed spherically 
symmetric mass clumps, and a Gaussian random-phase field. 

In both the scenarios we limit ourselves to large scale structures only, and do not attempt to 
include the lensing effects of smaller, more compact objects like individual galaxies and galaxy 
clusters. The computed correlation amplitude will not be affected by this exclusion: Theoretical 
calculations of weak lensing induced correlations effectively smooth out density fluctuations on 
scales below the correlation scale. This is seen clearly from equation (13) of Dolag &: Bartelmann 
(1997), where the cutoff is imposed by the filter function F(k,4>). 

On large scales the linear growth of mass density fluctuations is an adequate assumption; 
therefore at every cosmic epoch we scale the amplitude of fluctuations by a(t) = 1 + z\. 

6.2.1. Model #1: Clumps 

This model for the mass distribution is motivated by the observation that the large scale 
structure looks like a network of walls, voids and filaments, i.e. coherent structures spanning tens 
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of Megaparsecs. Therefore we represent the mass fluctuations by a collection of many spherically 
symmetric clumps. We have experimented with a number of different density profiles, including 
non-singular isothermal sphere, top hat, and Hernquist profile, p oc [r(l + r/ro) 3 ] -1 . Qualitatively 
the results described below are similar for all types of mass clump profiles provided the physical 
scale of the clumps is about the same. In what follows we use clumps described by a Newtonian 
potential, $(r) oc — e~ r / rc , where r is comoving distance from the clump center, and r c is its 
scale length. This particular shape was chosen because the potential falls off quickly so that the 
influence of the mass clumps does not extend over scales comparable to the Hubble length. The 
mass profile of a clump is 

p~ r / r c / f \ 

p( r)(X f ( 2 --. (21) 

The profile has a central mass peak (r < 2r c ) and is surrounded by a shallow region of the opposite 
density. The clumps are quite compact: for r c = 0.1-^ the density drops by a factor of ~ 50 
between r c /10 and r c , and by a factor of ~ 20 between r c and 3r c . The mass enclosed within r of 
each clump goes as e~ r ^ Tc r 2 r~ 1 , and asymptotically approaches zero. We use clumps of positive 
and negative central mass density, so even in a small volume the statistical expectation value for 
the net excess mass density is zero. This is in accord with the equations derived in Sections 2, 3 
and 4 which assume a flat space geometry. 

The clumps are randomly distributed in space with a typical interclump separation of about 
70/i~ 1 Mpc. The line of sight from the observer to the source is completely immersed in a sea 
of clumps. All clumps are identical, with r c = 0.1-^ = 300/i _1 Mpc which is larger than the 
interclump separation. Because of this the typical radius of the net structures formed by the 
superposition of clumps is about 0.05-^ = 150/i -1 Mpc, comparable to the size of the largest 
superclusters seen in the nearby Universe (Batuski et al. 1999). The rms fluctuations of mass 
in spheres of R = 300/i _1 Mpc is 0.012, comparable to, or somewhat higher than the standard 
ag = 1 normalized CDM value. The amplitude of fluctuations in this model grows slowly with 
decreasing scale; the fractional contribution to as is small. We will return to the discussion of the 
amplitude of mass fluctuations later, in Section 7. In the redshift range Z\ = 0.1 — 0.3, the location 
of the galaxies used in the cross-correlation studies, the average projected separation between 
two adjacent clumps is 25/i~ 1 Mpc, or about 2.5°, comparable to the angular scale of observed 
QSO-galaxy cross-correlations. 



6.2.2. Model #2: Gaussian random field 

Here we assume that the fluctuations are represented by a Gaussian random field. The 
form of the power spectrum was taken form Peacock (1997). To make a fair comparison 
with the clumps model we limit the power spectrum to a range of wavenumbers between 
how = 27r(2r c )~ 1 = 0.01/iMpc" 1 and khigh = 27r(0.5r c )~ 1 = 0.04/iMpc _1 . In this region the shape 
of the power spectrum is P(k) oc k. The phases of the various modes are uncorrelated. 
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7. Results of numerical integrations 

7.1. General considerations 

For every line of sight numerical integration of eq.(12) yields source magnification. Since we 
are interested in the cross-correlation between sources and the lensing mass, we first need to relate 
the source magnifications to the cross-correlation function. 

Williams and Irwin (1998) derived a relation, applicable in the weak lensing regime, between 
QSO-galaxy correlation function, ujqg(O), and the autocorrelation function of the foreground 
galaxies used in the analysis, ujgg{@) (their eq.[6]): ujqg{0) ~ (2r/6)(2.5a — 1)logg(9)i where r is 
the total optical depth of the lensing slab of matter traced by the APM galaxies, and b is the bias 
parameter. The factor (2r/6) is the proportionality coefficient between (M — 1), 'magnification 
excess' of a source, and normalized galaxy number excess, (cr — 1), both M and a referring to the 
same patch of the sky. In the weak lensing regime (M — 1) and (a — 1) are linearly related in 
both the standard and revised lensing formulations, but the coefficient is different. Thus for either 
formulation we can write 

ojQGiO) « ((M - 1)/(<t - l))(2.5a - l)u GG {6), (22) 

where the brackets indicate the typical value over many lines of sight. Without evaluating loqg{9) 
itself, eq.(22) allows us to derive an expression for the change in the predicted correlation function 
resulting from adding the second order term to eq.(4), while keeping the galaxy distribution, 
p(a\9) and the slope of source number counts, a the same: 

VQG,2nd(6) = / (M 2n d ~ 1) \ 
coqgM ) ^ ( M std ~ 1) / 

The quantity (Mmd — l)/{M s td — 1) for any line of sight has an interesting property of being 
independent of the amplitude of mass density fluctuations, while being sensitive to the shape of the 
fluctuations. This is readily seen from eq.(19) and (20): both versions of (M — 1) are proportional 
to $. Therefore the change in ujqg is the same regardless of the actual amplitude of fluctuations 
or the magnifications induced by the mass structures, so long as their shape remains the same. 

Although the actual values of magnifications are not important, the values presented below 
are comparable to what one expects from the smoothed large scale structure, i.e. <1%. 

All the integrations described below conserve flux. Flux is considered to be conserved if the 
ratio of the total amount of flux received by all observers distributed in a spherical shell around a 
source to the total flux emitted by that source is 1. For any given observer that ratio is just 1/M; 
the average 1/M over a 1000 or so observers in each set of integrations was found to be less than 
one standard deviation of the mean away from unity. 

We now describe the results for each of the two mass fluctuation models. 



(23) 
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7.2. Mass density fluctuations 

7.2.1. Model #1: Clumps 

Each point in the top panel of Figure 1 represents the magnification of a source computed 
using the standard and second order formulations: M st d and M2 n d, respectively. It is apparent 
that including the second order term has a small but non-negligible effect on the computed 
magnifications. The bottom panel, which shows the residuals AM = M 2n d — M st( i vs. M st d, 
demonstrates that M 2n d is, on the average, more extreme than M st d, in other words, the second 
order term increases the computed magnification of the source if M st d > 1, and decreases it if 
M st d < 1- This effect was already noted in Section 5. The fractional change is not very large; on 
the average AM/(M std - 1) « 0.2. 

Do these results alter the predicted amplitude of QSO-galaxy correlation function? Figure 2 
shows {Mmd — 1)/ {M st d — 1) as a function of the projected angular distance between the source and 
the nearest mass clump of either positive or negative central density. The only clumps considered 
here are in the redshift range 0.1 — 0.3, to simulate the observed situation. The angular scale of 
correlations includes the observational scale of <1°. The median of all the (M 2n d — l)/(M st d — 1) 
points is 1.1, i.e. the amplitude of the cross-correlation function is increased by 10% compared to 
the standard lensing equation. As the reader will remember from the Introduction, an increase by 
a factor of >5 is needed to reconcile observations with theory. 

The solid points in the plot represent sources (lines of sight) with M 2n d > 1-001, an arbitrary 
cutoff designed to separate out about half the magnification values. The median of the solid points 
is very similar to the one quoted above, i.e. 1.1, but the dispersion in the points is much smaller, 
and so it is more obvious that the points prefer to lie above the (M 2n d — l)/(Mst<2 ~~ 1) = 1 li ne - 
In a simple case where the source number counts truncate at a flux just below our detection limit 
the solid points will represent a flux limited sample of observed sources. 

Note that about half of all the points in this plot represent M st d < 1, i.e. demagnifications. In 
such cases AM values tend to be less than zero, therefore (M 2n d — l)/(M st d — 1) > 1, and lensing 
induced anti-correlations are also enhanced by the second order term. 

As was pointed out in Section 5 with the help of eq.(20), the presence of the second order term 
in eq.(4) may change the redshift of the optimal lenses compared to the standard formulation; it 
may 'move' the most effective lenses, for fixed source redshift, closer to the observer. The numerical 
integrations indicate such an effect. We have divided the nearby universe into four slices oriented 
perpendicular to the line of sight to the sources. The slices are: z\ = 0.05 — > 0.1, 0.1 — > 0.2, 
0.2 — > 0.3, and 0.3 — > 0.4. In each slice we compute the median value of {M 2n d — l)/(M st d — 1) for 
sources found within angle 9 of at least one clump. Only clumps of positive central density were 
used in this exercise to emulate the observations where QSOs are (anti) correlated with regions of 
galaxy excess. The results are plotted in Figure 3. Angle 9 was set to 0.5° (solid), 0.75° (dotted), 
and 1° (dashed), comparable to the angular scale of observed correlations. Mass concentrations 
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in all redshift slices are more correlated with the background sources than the standard lensing 
formalism would predict, however, the increase in the correlation amplitude is larger for z\ < 0.2 
slices than for the z\ > 0.2 slices, i.e. the inclusion of the second order term moves the redshift of 
the optimal (most correlated) mass fluctuations closer to the observer. 

How do the results change if we adopt a different value for the clump length scale, r c ? The 
effect of the second order term, evaluated with (M2 n d — l)/(Mstd — 1) declines rapidly as r c is 
reduced by a factor of 3 or larger. This is easy to understand; the effects of small scale density 
fluctuations on the net d^/dz cancel out as eq.(20) is integrated between the source and the 
observer. Only the effects of larger scale fluctuations, however small in amplitude, will remain. 

7.2.2. Model #2: Gaussian random field 

Figure 4 compares M st d with M^nd for the Gaussian random field mass fluctuation prescription, 
and is equivalent to Figure 1 of the clumps model. Though Gaussian random field does respond 
to the second order term, the change in magnification is very much smaller than in the clumps 
scenario, therefore there will be virtually no increase, over the standard lensing formulation, in 
the predicted amplitude of cross-correlations if a Gaussian random-phase field describes the mass 
density fluctuations on large length scales. 

Why does the inclusion of the second order term make a difference in the weak lensing model 
predictions when a clumps density model is used, while no effect is seen when the density field is 
Gaussian? As we have already mentioned the normalization of the power spectrum has no effect 
on (ujQG,2nd/^QG,std)- If normalization is increased the magnification labels on both the axis of 
Figure 4a will change, but the slope of the line defined by the points will not. Same is true for 
Figure la of the clumps model. The physical scale of density fluctuations is the same in both 
the models, so that too can be ruled out as the source of the difference in the results. The exact 
form of the clumps is also not a factor; we have tried a number of profile shapes, and all yielded 
very similar results. We are led to conclude that the difference is due to the shape of the density 
fluctuations in the clumps vs. Gaussian random field models, in particular it is due to the large 
scale coherent, i.e. non-Gaussian nature of the clumps. 

One can arrive at the same conclusion by considering the physical content of eqs.(ll) and (4). 
According to eq.(4) the second order term introduces some extra bending of the light rays towards 
positive mass inhomegeneities and away from negative ones. In a density field dominated by 
small-scale incoherent structures, like our Gaussian random field, the rays are bent many times in 
opposing directions as they travel from the source to the observer, thereby averaging to a straight 
line, a situation well described by the standard lensing approach. 
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8. Summary and Conclusions 

We have extended the standard lensing formulation to include the next order term (second 
order term in the geodesic equation) in the hope that it will make a significant difference in the 
predicted amplitude of weak lensing induced QSO-galaxy correlations, which is currently at least 
a factor of 5 smaller than observed. Depending on the type of the mass density fluctuations 
populating our model universe, we can obtain up to a ~ 10% increase in the amplitude of the 
correlations, i.e. not enough to reconcile theory with existing observations. 

Even though the increase is small, qualitative aspects of the updated weak lensing formulation 
bring us closer to the properties of the observed correlations. The revised source magnifications 
correlate with, but are more extreme than the magnifications computed using the standard 
formulation, i.e. typically (M2 n d — l)/(Mst<2 — 1) is greater than 1, which increases the amplitude 
of predicted lensing-induced correlations. Furthermore, the redshift of the optimal lenses for 
cross-correlations with sources at a given redshift is closer to the observer than what the standard 
formulation would predict. This is consistent with observational results. It is important to note 
that the second order term affects only the predictions of the lensing theory which deal with 
weak statistical lensing of point sources. As we have shown using eq.(17) and (20) and numerical 
integrations, image elongation, i.e. weak shear lensing, and strong lensing theory are not altered. 

Potentially the most important result of the present work is that not all types of cosmic 
density fluctuations are sensitive to the presence of the next higher order term in the lensing 
equation. Only large-scale, coherent, i.e. non-Gaussian structures show an increase in ujqg with 
the revised formulation. This can be intuitively understood in terms of eq.(ll) and (4) (see end 
of previous Section), and quantified with eq.(20), which states that the increase in magnification 
induced by the second order term is proportional to a weighted integral of d<&/dz over the line of 
sight. The integral will tend to be smaller when the phases of the modes comprising the density 
and the potential field are uncorrelated, compared to a density fluctuation model dominated by 
coherent structures. Note that the increase in the predicted cross-correlation, (uJQG,2nd/^>QG,std) 
does not depend on the amplitude of the mass fluctuations, only on their form. 

The work presented here is preliminary, however, if the avenue we have taken proves to be 
correct, i.e. if second order, and possibly higher order terms in the geodesic equation of motion of 
photons are important for weak lensing induced QSO-galaxy correlations, then these correlations 
will give us the means to study the type of mass structures present on the largest scales, in 
particular to test the Gaussianity of the large scale fluctuations and thus shed light on the physical 
processes in the early Universe that gave rise to the present day structures. 

I am deeply indebted to Prasenjit Saha for numerous invaluable conversations on the subject 
and many helpful suggestions. I am grateful to the anonymous referee for further elucidating the 
physical meaning of the second order term. 
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Fig. 1. — Magnification of sources obtained using the standard lensing formalism, M stc i, and 
magnification computed when the second order term in the geodesic equation is included, Mmd- 
The mass fluctuations in this model universe are build up with a collection of randomly distributed 
spherical clumps of positive and negative central density. See Section 6 for details. The residuals, 
AM, plotted in the bottom panel, are positively correlated with M stt ±, i.e. inclusion of the second 
order term enhances (de)magnifications. 
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Fig. 2. — The factor by which cross-correlations are increased is plotted against the projected 
distance to the nearest mass clump in the z\ = 0.1 — 0.3 redshift range. All the sources are at 
z s = 1.5, and mass density fluctuations are represented by a collection of spherically symmetric 
mass clumps (see Section 6 for more details.) The median increase in the amplitude of uiqg is 
1.1. The solid points are those with M2 n d > 1-001 and have the same median. Correlations are 
increased on all angular scales. 
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Fig. 3. — Change in the amplitude of predicted cross-correlation function between z s = 1.5 sources 
and mass distribution in four low-redshift slices, z\ = 0.05 — > 0.1, 0.1 — > 0.2, 0.2 — > 0.3, and 
0.3 — ► 0.4. In each of the four redshift slices the median value of ((M2 n d — l)/(Mstd — !))> equal 
to {ojQG,2nd/^>QG,std)i was computed for sources lying within of a mass clump (positive central 
mass density clumps were used here). Angle 9 was chosen to be comparable to the angular scale of 
observations, 6 = 0.5° (solid), 0.75° (dotted), and 1.0° (dashed). (Note that the number of lines of 
sight contributing to the twelve points of the plot varies significantly: the = 0.5°, z\ = 0.05 — ► 0.1 
point has 20 lines of sight, while = 1.0°, 0.3 — > 0.4 point has ~ 1700.) The inclusion of the second 
order term in the geodesic equation preferentially increases the amplitude of correlations for more 
nearby lenses, thus 'moving' the optimal lenses closer to the observer, compared to the standard 
optimal lensing distance. 
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Fig. 4. — Similar to Figure 1, except the mass fluctuations here are described by a Gaussian 
random-phase field, with the power spectrum, P(k) oc k, and limited to the range of wavenumbers 
between k tow = 2vr(2r c )" 1 = 0.01/iMpc" 1 and k high = 2vr(0.5r c )" 1 = 0.04/iMpc" 1 . The residuals, 
AM are much smaller than in Figure 1. 



